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Abstract The soft wall model in holographic QCD has Regge trajectories but wrong oper¬ 
ator product expansion (OPE) for the two-point vectorial QCD Green function. We modify 
the dilaton potential to comply with the OPE. We study also the axial two-point function us¬ 
ing the same modified dilaton field and an additional scalar field to address chiral symmetry 
breaking. OPE is recovered adding a boundary term and low energy chiral parameters, 
and Lio, are well described analytically by the model in terms of Regge spacing and QCD 
condensates. The model nicely supports and extends previous theoretical analyses advocat¬ 
ing Digamma function to study QCD two-point functions in different momentum regions. 


1 Introduction 

QCD Green functions describe three ranges of energies (i) deep euclidean, where pertur¬ 
bative QCD and OPE expansion methods can be applied, (ii) an intermediate Minkowski 
region, where resonances are described by Regge trajectories and (Hi) the strong interacting, 
low-energy region, described by chiral perturbation theory (^PT) [1]. 

Regge trajectories, i.e. a linear growth of the square of the resonance masses like My [n) = 
On GeV^ for the vectors, were conjectured long ago, before QCD, and very well verified 
phenomenologically [2]. This relation, thought to be valid in QCD for large number of 
colours, Nc, has been proven by’t Hooft only in 2 dimensions [3]. Also other arguments 
have been proposed for the validity of Regge theory aside from the experimental evidence: 
lattice, flux tube, confinement [4-6]: we insist on these theoretical motivations of Regge 
trajectories through this paper since we believe that QCD dynamics favours this physical 
picture. 

Low energy QCD properties, like the chiral symmetry breaking (^SB) parameters, 
and the Gasser-Leutwyler coefficients L, ’s of the chiral Lagrangian have been studied 

in the framework of holographic models of QCD, based on the AdS/CFT correspondence 
[7-9]. Several five dimensional set-ups have been proposed. The basic feature of Hard-Wall 
(HW) models is to simulate confinement by cutting drastically the extradimension of the 
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AdSs in the IR, producing an infinite spectrum of Kaluza-Klein (KK) states to be identified 
with vector resonances of increasing masses. ;^SB in the axial sector is triggered either by 
a scalar field in the bulk [10, 11] or by appropriate boundary conditions [12]. These models 
were actually anticipated by a deconstruction model [13]. Also the Sakai-Sugimoto model 

[14] , a string set- up, shares this feature of an infinite spectrum of KK states. A recent review 
of all these models and their relation with the light-front holographic QCD approach is Ref 

[15] . 

All these models have a good description of the deep euclidean region; for instance the 
correct two point function recover the partonic log, and a pretty well description of low 
energy QCD, obtaining the chiral parameters and the Gasser-Leutwyler coefficients L/’s 
which are close to the physical values. In fact, it was observed that these HW models [10-12] 
have the same two point vectorial Green function as the one obtained by Migdal long ago 
who proposed an ad hoc prescription to perform the analytic continuation of the perturbative 
deep Euclidean QCD result to the Minkowski region [16, 17]. 

In order to reproduce resonances masses with a Regge spacing, one can consider a 5D 
model with AdS metric and an additional field, the dilaton [18]. However, in this model it 
can be shown that the partonic log of the two point vectorial Green function 



( 1 ) 


with the Euclidean momentum = —q^, receives '/e^-corrections differently from what 
holds in perturbative QCD Operator Product Expansion (OPE)* [19-22]. Nevertheless, we 
find extremely interesting to have a model where the Regge region would be analytically 
related to the deep Euclidean region even though with the wrong OPE. This will be exactly 
our starting point: OPE tells us the correct Green functions in the deep Euclidean in terms 
of gluon and quarks condensates. Is it possible to modify the dilaton profile, (j){z) (j) (z) + 

5(j){z), such to comply with QCD requirements in the intermediate (Regge) region and UV 
(OPE) region? Once this question has been answered in the affirmative, we are naturally led 
to add ress another important problem, that does not have a satisfying answer in this "linear 
confinement" holographic approach [18]: how to implement the effects of chiral symmetry 
breaking into our model. If we look at the expansion of both vector-vector and axial-axial 
correlator TIvAiQ^) fo'" large Euclidean momentum, QCD predicts 



(2) 


with cy = 1 and ca = — "/v. We should then be able to take into account the difference 
appearing at order '/e'^ in the OPE expansion of the vector and axial correlators. On the other 
hand, in the low and intermediate region in Minkowskian momentum, chiral symmetry is 
broken at the level of the mass spectra of resonances, with the appearance of the pion as the 
corresponding pseudo-Goldston boson of ^SB. If we were able to deal with these aspects 
in a holographic model, we would have, in principle, closed the circle, since we would have 
a description of all three energy ranges of QCD. In fact, implementing the correct OPE and 
mass spectrum for the axial sector, we will propose our post-diction of and Cio (one of 
the the chiral coefficients) in terms of our input parameters, i.e. the Regge spacing 

and the QCD condensates. Our model, beside being a novel proposal for holographic QCD, 


*Due to its definition, Tly used in the holographic models [10] differs by a factor '/i from the fly used in 
[1,4]. 
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links naturally to previous 4D QCD work where a phenomenological matching between low 
energy and OPE was realised through Regge theory (or Veneziano model) [1,4, 23-27]. 

We have organized this paper in the following way. In Section 2 we summarize the 
relevant properties of the existing Soft Wall model (SW) and focus on the vectorial two-point 
function. The issue of ^SB is dealt with in Section 3. In section 4, the analytic continuation 
at low energy is described and the predictions for and Lio are given as well as some 
consequences on duality violations. Conclusions are given in section 5. 


2 Assumptions and aims for the model 

2.1 Preliminaries 

The Soft-Wall model is a five dimensional model where the additional coordinate, z, has 
the range 0 < z < °° and background fields consist in a five-dimensional AdS metric and a 
dilaton field 4’(z). The AdS metric is written as 

gMN^^ =-^ {r|^lvdx^^''-dz^) , (3) 

where ri^v = Diag(l, — 1, — 1, —1), the Greek indices It, V = (0,1,2,3) referring to the usual 
4-dimensions, and the capital Latin ones M,N = (0, l,2,3,z) to the 5 dimensions. It was 
shown, in [18], that with the choice of a quadratic profile of the dilaton field 

^{z) = K-V, (4) 

the spectrum of vector resonances followed a Regge trajectory. Vector meson resonances are 
obtained as the modes associated to a five-dimensional gauge field Vm in the external metric 
and dilaton background, with 5D action 

55 = |d4xjr”dzVJe-‘*’W Tr[FM«F^5] , (5) 

with g = I detgMwl, the field strength F^w = dt/V m — m — i [Vm, Va^] and gj = is 

the 5D coupling constant where Nc is the number of colors of QCD. The trace is understood 
as the sum over the S\J{Nf) flavour indices ofY^ = if Nf = 2, t‘’ = ^ ji with being 

Pauli matrices. 

We shall work in the axial gauge = 0. The AdS/CFT correspondence prescribes that 
the boundary value of the 5D gauge field has to be identified with the classical source 
coupled to the the 4-dimensional vectorial current 7y ='YYii h 

= . ( 6 ) 

The corresponding equation of motion for the gauge field, derived from the Lagrangian 
(5), is more easily written in terms of the 4-dimensional Fourier transform fv{—‘f',z) of the 
field Vm, 

fv{-q^,z) = j d'^xe^“^'^Yfifi{x,z) , (7) 

where is the Fourier transform of the source v^j. Using gauge invariance, one may assume 
the Vfi to be transverse and one obtains 


(9| fv + [In wo (z)] d^fv+q^fv=0, 


(8) 



4 


with, 

wo(z) =-=-■ (9) 

z z 

Notice that, in order to satisfy the the boundary condition (6), we have to impose 

/v(V,0) = l. (10) 

It is well known that the presence of (an infinite set of) solutions of the 5D equation 
of motion (8), satisfying a normalization condition in the extra-dimension z, corresponds 
to (an infinite set of) resonances in 4D. Indeed, the corresponding eigenvalue equation for 
wave functions (^„(z) of the normalizable modes can be obtained from (8) by requiring 
= Mv{n)^, which gives (^„(z) = fv{—Mv{n)^,z) and the equation 

- d^(l)„+Mv{nf(l)„=0, ( 11 ) 

with the vanishing boundary conditions 

(t>n{0)=0 and (12) 

The authors of Ref. [18] showed that, by doing the substitution 

(?„(z) = \/z%{z'), (13) 

one transforms the equation (11) into a Schrodinger equation for 'f^(z) 

+ (14) 

which is exactly solvable in terms of generalized Laguerre polynomials L„ . The solutions 

(^„(z) of the original equation (11) are then 

(?„(z) = K^Z^J , (15) 

V n+ 1 

and the required normalization conditions are given for each representations by 

[ dzwo{z)^m{z)^n{z) = 5,„„ and [ dz'I',„{z)%{z) = 5,„„ . (16) 

Jo Jo 

The most important result, which characterizes the Soft-Wall model and makes it a pos¬ 
sible improvement, with respect to the HW models, concerns the resulting mass spectrum 
which was found in [18] to be 

Mvinf = 4K^{n+l) (17) 

with w = 0,1,_Eq. (17) shows that the infinite set of 4D resonances corresponding to the 

5D eigenfunctions follows a Regge trajectory. 

Despite this very nice approach to understand the Regge trajectory, the SW model does 
not reproduce the correct OPE of the vector correlator. In fact, the large-Q^ expansion of 
two-point function coming from the solution of (8) has a non vanishing cofficient for the 
'/q- term, corresponding a non vanishing dimension-two condensate, in contrast to what 
predicted in QCD (a fact that was already pointed out in [19-21]). Also the values of higher 
dimension condensates turn out to be different from those obtained in QCD. 

The first purpose of this paper will be build a deformation of the dilaton field such that 
the OPE of the vector correlator corresponds to the real one. 
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2.2 The vector two-point function 


One defines the vector-vector correlator Tlv , 


ijd'^xe “>'^(^Jvii{x)Jyv{ 0)'^ = 5‘‘'’{qnqv-q'^rinv)nv{-q^) . (18) 

where q is the 4D momentum and the current Jy ^ is the one defined before equation (6). 

The function fv{—q^,z) satisfying eq. (8), (10) and vanishing at z ^ the so-called the 
bulk-to-boundary propagator, plays a pivotal role in holographic models, since its knowledge 
is fundamental in evaluating Green functions of the 4D theory. Once fy has been found, it 
can be substituted in the quadratic part of the 5D action (5), which is the one relevant for the 
evaluation of the two-point function. The resulting 5D expression reduces to a 4D boundary 
term, quadratic in the 4D source , from which one gets the two-point function as the limit 
(with the Euclidean momentum = —q^) 


Q^ny{Q^) = ^ limwo(z)/y(g^,z) dJy{Q^,z) ■ 


(19) 


Our main concern is that the two-point function Tly coming from the model built in 5D 
has the following properties: 


(i) The progression of its poles in the Minkowski region (the resonances) follows a 
Regge trajectory. 

In the Large-Nc limit of QCD, the two-point function Fly could be written as a sum 
over an infinite set of stable vector resonances [3 1 , 32] 


nv{Q^) 


y Fy{nf 
^oQ'^+My(n)^ ’ 


( 20 ) 


where Fy{n) are named decay constants and the My{ri) are the masses associated to 
the resonances of the vectorial channel: p, p', p". We assume that these resonances 
follow, in a first approximation, a Regge progression [2], 

My{nf ~ an, (21) 

«—^■oo 

where the integer n is the radial excitation number and a is related to the confining 
string tension as explained in [18] and we can evaluate a 0.90 GeV^ from [33] as 
illustrated in figure 1 . As pointed out in [34] by their authors, one has to organize the 
rho resonances according to their radial and their angular-momentum progressions in 
order to select the right Regge trajectory. It means that we should not consider all set 
of resonances in figure (la) but only the ones in figure (lb) corresponding to the 
masses squared progression My{n)^ ~ 1.43(13)n GeV^ or ct w 1.43(13) GeV^ One 
assumes here that the Regge behaviour is more and more valid the larger n is, which, 
of course, implies that the description of the lower resonances contributions would be 
less precise. From now on this spectrum will be called the physical spectrum. 

(ii) Fly has the correct OPE. 

One of the very well-known properties of this two point function is its OPE [35], 


nv{Q^) 


1 Nr 


22^00 2 12;r2 


In 



1 


1 


{^6)y 
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^We thank Pere Masjuan to have pointed out this observation. 


( 22 ) 
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Fig. 1 The progression of the squared masses of the p resonances according to their label n. The straight line 
represents the fitted relation. The (a) plot follows [18] for My[nf- ~ 0.90« GeV^. The (b) plot represents the 
fit obtained in [34] with My{n)^ ~ 1.43(13)« GeV^. 


where in the large-Nc limit the coefficients of the OPE are given by 

here a, stands for the strong coupling constant, (G^) for the gluon condensate, and 
{yy) for the quark condensate. 

The requirement that the two-point function obtained from the SW model (5) with a 
suitable deformation of the dilaton profile, has the correct OPE (22) and (23), is the starting 
point of our analysis. Proposals to modify the original SW model based on the assumption 
of different profiles for the dilaton and for the vacuum expectation of a scalar field in the 
bulk were done in [28, 29]. 


2.3 Vectorial OPE from a modified dilatonic background 

We shall assume that the effects of the OPE on the vector current two-point functions in QCD 
can be encoded in a new profile for the dilaton field of the SW model. More explicitly, we 
assume that OPE is related to the behaviour of the dilaton profile around the UV boundary, 
z = 0 of the 5D, while keeping the metric to have the AdS form (3). Moreover, as we want 
to keep the Regge behaviour induced by the dilaton profile 4’(z), we are led to modify the 
original quadratic profile of the SW dilaton by adding new terms which we collect in a 
function B, 

4>(z)^<P(z)+5(z). (24) 

We assume that the function B can be represented for all z by a polynomial of degree 
2K (with no constant term), 

k=l 

Even for a polynomial form of the dilaton profile as given in (1), the task of solving the 
corresponding equation of motion of the vector field is a non trivial one. In the following, 
we show how the ad hoc peculiar form for 5(z) allows us to solve the equation of motion 
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analytically, through an iterative method, and obtain the corresponding expression for the 
vector two-point function. To simplify the writing we define, 

wo(z) = w[z) ■ (26) 

For counting the different orders of perturbation of the dilaton, we introduce an artificial 
control parameter Q such as 

w(z) = . (27) 

We shall expand all relevant expressions in (formal) power series in 0, using this parameter 
to translate the task of solving the differential equation for the vector field in the deformed 
dilatonic background in a more tractable set of iterative equations. 

At the end of the calculation, we shall match the asymptotic expansion of the two-point 
function to the OPE expansion of QCD. Doing this the dependence of the solution on the 
parameter 0 disappears. We shall see that each term in the expression of B{z) drives to 
a term in the OPE of the two point function and since we want to match the OPE only 
up to '/s'* term, we shall need only a polynomial of degree 6, i.e. to take A" = 3 in (25). The 
matching with the QCD OPE uniquely determines the coefficients of the dilaton profile. 

Let notice that, of course, other forms for B{z) could give the same OPE but different 
exact solutions of the equation of motion to be interpreted as analytic continuation of the 
OPE. Our expression for B is the simplest way to satisfy our assumptions ^ . 

In the next Section 2.4 we explain the technical details of the iterative method we derived 
to obtain the vector two-point function for the modified dilaton background. The complete 
expression for the two-point function is in Appendix B. 


2.4 Iterative construction of the vector two-point function for the modified dilaton 
background 

As we said before, if we want to reproduce only up to the term, the calculation has to 
be done up to order 0^. 

Using (27), the equation of motion for the vectorial field V associated to wq when ex¬ 
pressed in its Fourier transform fy becomes, 

fv + dz [Inw(z)] dzfv-Q^ fv = 0 , (28) 

and can be decomposed into a hierarchical system of differential equations order by order 
in 0. 

We can expand fy in power of 0, 

fv = + Qfv^ + 0 Vf + e + (29) 

with the boundary conditions 

/W(e2,0) = 5o,„ and /W(g2,c»)=0. (30) 

^One could wonder if a complete knowledge of the full OPE expansion would allow us to write B(z) as an 
infinite power series. Even in such an ideal case, that series would probably be only asymptotic and additional 
dynamical assumptions on how B extends at any value of z would be unavoidable, in the same manner that 
we have to deal with renormalons prescriptions [30]. 
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Then, the associated two-point function (19) has an expansion in 0 too, 

Hv = -h -p ff{e^) , (3i) 

where 

Q^n^\Q^) = 4 limwo(z)4°^(e2,z)<9,4'^\e2^z), (32) 

thanks to the boundary conditions (30). 

To simplify a little the writing, let us define the differential operator 33 such as 

Dtp = 5^9-b <9;. [In wo(z)]<9z(p-gV (33) 

then the equation 'Sfy'^ = 0 corresponds to the unperturbed equation of motion (8). Or¬ 
ganising order by order in 0 the differential equations, one obtains the recursive system. 


= E ‘ b2in-k) d, /f , (34) 

k=0 

or explicitly, for the first three orders, 

23/^ = 0 (35) 

Dfl}^=zb2d,4°^ (36) 

D/f = b4 + z b 2 <9,/^^ (37) 

be b4 -b Z ^2 <9z/y^'. (38) 

To solve the system, we can use the Green function method. Indeed, from the first order 
correction (36) on, one has the same differential equation with different source terms. To 
solve any equation with a source term S like 

T)fv = S{Q\z), (39) 

one first solves the equation for the Green function Gv{Q^\x,z), 

T)Gv{Q^',x,z) = 5{x-z), (40) 


with the boundaries conditions, 

Gv{Q^Az) = Gv{Q^;^,z) = 0. 

Then, the solution of (39) is given by 

fv{Q^,z)= [ dxwo(x) Gv(Q^-,z,x) S(Q^,x) . 

Jo 


(41) 


(42) 


Therefore, the differential system has for solution the recursive expression (for n > 0), 

'n-l 


fv\Q^,z)= f dxwo{x) Gv{Q^-,x,z) 
Jo 


^ b2(n-k) dxfv\Q^,x) 


k=0 


(43) 


and using (32), the solution for the vector two-point function is the convoluted expression, 
b2(n^k) [ ^2{n-k-i) dxfl 2 ^\Q^,x) , (44) 

^5 to Jo 
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where we used the relation, 

i°\Q\x) = -liInwo(^)/^(e^^)^^Gy(e^x,^) . (45) 


In order to solve this system recursively, we have to know first fy \ which is nothing 
else that the well-known solution of the unperturbed equation of motion 'Sfy^ = 0 i.e. the 
original bulk-to-boundary propagator of the SW model. It could be expressed in several 
ways, for instance, as a series of the eigenfunctions, which in the SW model are given in 
terms of Laguerre polynomials (15) [18] . 

However, in order to perform the calculation, we found more efficient to use the integral 
representation for flp given in [36], i.e. 


AO) 

JV 




du M4k- 


4- 


exp 


1 — M 


2 2 
K^Z 


(46) 


moreover, we constructed a similar representation also for the Green function (cf. Appendix A), 


Gv{Q^-,x,y) = - y ^ df 


e- _ 1 

f ii? 3 
l-t 


■exp 


l-t 


K^{x^+y^) 


Ii iK^^xy 


l-t 


Leaving the details of the calculation to the Appendix B, the final result for Tly'^ is 




1-f 




\4(c2 / \ 4 k ^ ) 


4(c2 / 


(47) 


(48) 


in) 

In the same way, one can express order by order the analytic expression for Tly which takes 
the general form 

Q^n^\Q^) = (£) (^) , ( 49 ) 

where, respectively, !^k,n are polynomials and is the k* derivative of the Digamma iff 
function defined as the logarithmic derivative of the Euler’s F function. The coefficients of 
„ depend only on and the coefficients bi^ of the dilaton (25). The explicit expressions 
for and are given respectively in eq. (B.19) and (B.20) in Appendix B. 

Our method provides a framework where we can get an explicit solution in terms of hy¬ 
pergeometric functions (see eqs. (B.15) and (B.16) in Appendix B), which then reduce to a 
combination of Digamma functions, and, importantly, it avoids the appearance of divergent 
series and divergent integrals, as it can be seen for instance, in the detailed calculation of 
riy'^ presented in Appendix B. Moreover, as shown in Appendix E, the use of such inte¬ 
gral representations allows us to organize the calculation of each order correction in a very 
systematic way. 

Let us comment here that the function. 




2 k -2 f Q^\ 
gj \4 k^) 


YE + Xj/ 



(50) 


can be associated to the ad hoc models developed in [24, 27, 37] with the same Regge 
progression than us. From the Regge spectrum point of view itself, is also similar 
to expressions developed in [38] (see references therein) obtained by resummation over 


resonances. 
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2.5 Interpretation of the solution 

Given the complete analytic solution of Tly, 


n=0 


(51) 


we can ask ourself if (i) the poles ofTly in the Minkowski region still follow a Regge trajec¬ 
tory and (ii) a correct OPE is obtained. We find that: 

In) 

(i) Since the part of each Tlf containing the Digamma function and its derivatives have 
poles only at all negative integers (—n), 


-d 


2 


4(c2 


= -«, 


(52) 


at 0-th order, only simple poles appear in Uy^'^, and so the usual Regge spectrum (notic¬ 
ing also that the residue in n = 0 is nought). 


M(n)^ = 4lf^n = an , 


(53) 


is reproduced, schematically presented in figure 2. However, in eq. (51), poles of higher 
order appear and a more detailed analysis has to be done to quantify possible departures 
from the original Regge spectrum. 


0 12 3 4 

n 


Fig. 2 Vectorial Regge spectrum 


(ii) In order to reproduce the correct OPE up to the order i/o‘, one can perform the asymp¬ 
totic expansion for large of ny(Q^) from the analytic expressions for the in 
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Appendix B, 


nv{Q^) 


' 4ic^e 


=2^2-"V a 




^mr 


^ogi 

4 

L V-- ■ ^4 , ^o;j g6 


+ ^ [- (2K:2 + eZ72)e^^4 +40^^76)] 


and match (54) term by term to the expected OPE 


nv(Q^) 


Nc 

247Z^ 






712 + ^^ 2 ) 712 + (^ 4 ) TjT + {^6)v Tjz 


Q 


G^ 


G® 


(54) 


(55) 


with the values of the condensates given in (23). Notice that Ay is automatically fixed 

_ '^E 

to Iks t ft! 1 GeV, which is the natural cut-off for our model. It implies also that 

'b2 = -2K^ 

< h4 = yi(ff4) (56) 

b6 = ^g5{^6)v 

which completely fixes w(z) and hence the dilaton profile B{z) in (25). In (56) we have 
removed the 6 dependence since if one uses these expressions in B{\fdz), the artifi¬ 
cial parameter 0 disappears naturally. From now, we will no more make explicit the 
dependence in 0 to simplify expressions for the reader except when this dependence 
is useful for the comprehension. 

Let us notice here, that as we anticipated in Sect. 2.3, (56) shows that there is one-by-one 
correspondence between a z^^ term in the function B and a term at '/q“ in the OPE. 


2.6 Analysis of the vectorial spectrum 

From the complete and analytic expression for Tlv, one has shown that the poles of the 
spectrum follow a Regge behaviour, with masses squared M{n)^ = The Large-Nc 

representation of Fly in (20) implies that the residues associated to the masses poles define 
the so-called Large-Nc decay constants 

Fy [nf = Res [Uy (Q^) , -4K^n] , (57) 

with Fy{0) = 0. They are obtained from the analytic expressions of each JT^^^ given in 
Appendix B. Despite this canonical description of the vectorial correlator, the Fy{n)^ cannot 
fully reconstruct the analytic expression of Fly. For example, let us focus on the order 0^ 
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associated to the function Tl, 


( 2 ) 


V 


. As expressed in (B.19), 




b4 / 4if^ 

V 2 ^ 


+ 


- 2 - 


^2 


l6K^gl 


- 1+2 


4k:2 


5 + 6 


4k'2 


+ 2 


_^y 

4k^ ) 


V 4(c2 


V‘ 


\ 4K'2 


+ 


1+3 

4K'2 


4(c2 

2 


^ ' 4k-2 


I 4(c2 


(58) 


and the corresponding residue progression is 

Res (g^),= ^(1 -6n) , 
Ss 


(59) 


that manifestly does not contain the part proportional to in (58). Then, whatever sum¬ 
mation over resonances procedure is used, one cannot recover the complete function. More 
generally, the analytic expression (49) contains derivatives of the Digamma function, which, 
due to their expression 





n=0 



(60) 


cannot be directly written as series of single poles. Physically, the previous expression of 
Digamma derivatives corresponds to the combination of A:+ 1 propagators, this suggests that 
the analytic expression of the vectorial correlator is a combination of several propagators. 
If one looks at this problem from a "perturbation theory" point of view, it could be possible 
to re-encode the contributions from several propagators into only one propagator but with a 
modified mass. In our case, it means that the higher order corrections in 0, or equivalently in 
the OPE, could be expressed as subleading contribution to the Regge behaviour, as in Ref. 
[39], for more details see Appendix D. 


3 Axial two point function 

In the previous section, we have presented a SW model with a deformed dilaton profile 
which describes Regge theory and allows us to obtain the correct OPE of the vectorial two- 
point function. The natural question is now if it is possible to do a similar construction for the 
axial correlator. The problem to recover the right axial spectrum and the right axial OPE is 
more complicated if we assume, as usual, that the coupling of the 5D vector and axial vector 
gauge field to the metric and the dilaton is the same. Actually, under these assumptions, the 
form of the dilaton profile is already fixed by the requirements of a correct OPE of the vector 
two-point function, with the coefficients identified to the vectorial OPE coefficients (56). 
So we need new ingredients to reproduce the QCD patterns of chiral symmetry breaking, as 
we already mentioned at the beginning of Sect. 2.1. 

Let first enlarge the gauge symmetry of the 5D Lagrangian to the whole group of chiral 
transformations SU(A/)l x SU(A/)«, introducing gauge fields, denoted by and 

each transforming only under one of the two SU(A/) groups. Eurthermore, as in 
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[10, 11], we introduce also a 5D scalar field X(jc,z) transforming under the bi-fundamental 
of the chiral group so that the gauge invariant action is taken as 

S5 = l Jd^x 1^ Tr|g"'^(DMX)^ (D^X) |xp 

- ^mNgRS ) I , (61) 

where = <9"X - iL"X -f iXM", dV- /[L", L^] and analogous 

expression for F^^. 

The action, can be rewritten in terms of vector and axial fields, V = '/ 2 (L-|-]R.) and 
A = V 2 (L — M), and as shown in [10, 11], taking 

x« = ^r, (62) 

the equation of motion for the Fourier transform over the 4D space of the axial field A, 
/a(— becomes 

d^fA + d,[\nw{z)]djA-Q^fA = g5(^^^ fA, (63) 

while the equation of motion for the vector field remains unchanged. 

In this approach, chiral symmetry is broken by the 5D scalar field, and in particular 
by a non vanishing profile v(z). The form of v(z) near the origin z ~ 0 is dictated in the 
AdS/CFT correspondence by asking the field X to be dual to the bilinear quark qq operator, 
whose non vanishing VEV is responsible for spontaneous ;^SB in QCD. In the following, 
we shall assume for the dilaton profile B(z) the one obtained from the OPE coefficients of 
the vectorial fields and then use a suitable form of the profile v(z) to encode the properties 
of the axial sector: 

- The axial spectrum contains a pion pole at = 0 and a Regge spectrum starting at 
q^ = with the same spacing than the vectorial spectrum. 

- The axial OPE has the following expression 

. 2^ (i)^^ 

and 

(^6^ = -y(^6)v . (65) 

with the same definitions than in (23). 


3.1 Realization of the axial spectrum 


In order to satisfy the axial spectral properties, we make the following first ansatz for the 
contribution due to the scalar profile in (63): 


v(z) 

z 


ft + jS*z5(z). 


( 66 ) 




14 


0 12 3 4 

n 


Fig. 3 Axial Regge spectrum 


Let focus first on the Regge progression of the axial spectrum. We notice that, phe¬ 
nomenologically, taking the first axial pole at 2Mp ~ 2a is quite an acceptable 

approximation for the axial spectrum. Then our prescription to obtain the axial spectrum 
from the vectorial one is to perform a "shift" over the vectorial spectrum like Q^+4k^, 
as suggested from the comparison of the vectorial and axial mass spectra (see figures 2 and 
3). This fixes directly in (63) 

glPo = 4K^ = a . (67) 

The full axial spectrum is not yet realised, indeed, since the axial spectrum is a shifted 
version of the vectorial spectrum, it does not contain yet any pion pole and moreover on the 
Euclidean side, this shift implies a modifcation of the OPE as 


2gj 


4K^e \ 2k^ 1 




( 68 ) 


where absence of dimension two operator in the axial OPE is violated. 

The introduction of the Dirac delta function in (66) has the nice properties to cure this 
two problems at the same time. This term in (63) generates only one exact contribution 
—^ (for any Q^) such that the pion pole appears naturally and modifying the OPE as 




4if^e 


-Je 


G" 


+ 


2k^ 



(69) 


then by taking 

2ic^ 

P* = —, (70) 

^5 

the axial OPE properties is satisfied again. The existence of a pion wave function associated 
to this mechanism and its use in the evaluation of low energy chiral parameters will be 
exposed in details in a forthcoming paper [40]. 
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3.2 Realization of the axial OPE 


The generation of the rest of the axial OPE terms derives exactly from the same procedure 
used in the vectorial section provided that we add two other terms in the expression ( 66 ) 
such as 

(t^) =A)+jS*z5(z)+j62z2 + /34z^ (71) 

With slight modifications, the iterative method we used to obtain the corrections to the 
vector two-point function, can be applied to the axial case too. All the details can be found 
in Appendix C. Then one has by identification: = Ay, the coefficients P 2 and 

P 4 are fixed by matching with QCD axial OPE (64) ,(65) : 

6 jA 

(72) 



The polynomial part of the scalar profile v(z) is 


v(z) = zVjSo + + /34Z^ 

2k 3k^ 3 / 61k^ 

z-^o g5 ^ 2g5 ^ V 48g5 


SgsK- 

4 


(^ 4 ) + 


45g5 

112k: 



(73) 


where the first two terms, which are the leading terms near z = 0 , are exactly the ones of 
a scalar field dual to the bilinear quark qq operator, i.e. the one required in Ref. [28] and 
compatible with the approach in [41]. 

The additional distributional term, proportional to j 6 * in (71) can also be interpreted as a 
boundary term, such that in the 5D action it would correspond to a peculiar 4D source term 
in the generating functional, 

P*j(fxj dzz5(z)w{z) A^{x,z)^ = P* J (fxA^{x,Q) A^{x,0) . (74) 


The appearance of the Dirac 5 function term, dictated by the request of correct reso¬ 
nance mass spectrum and OPE for the axial two-point function, makes our model definitely 
different from the more usual approaches adopted for instance in [28] and [42]. 


4 Analytic continuation in the chiral sector: the left-right correlator 

In previous sections we have explicitly shown how to implement the constraints on an axial 
and vectorial two-point functions coming from two different regions in the 17 ^—plane: the 
deep Euclidean region where we reproduce the matching with the OPE of QCD, and the 
Minkowski region where the two-point function have poles following Regge trajectories. 
Having built explicit expressions for Fly and FIa valid on the whole complex plane we now 
turn to the analysis of their prediction for chiral quantities defined at low Q^. 
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4.1 The left-right correlator spectrum 

Since we are now interested to the chiral sector, i.e. the low expansions, it is more perti¬ 
nent to consider the TI/jj correlator. 


nLR(Q^)=^-{nviQ^)-nA{Q^)) . 


(75) 


which is an order parameter of the chiral symmetry breaking mechanism in QCD. 


4.2 Predictions for the chiral constants 

The low limit can be obtained by the analytic continuation of our expressions for the 
axial and vectorial two-point function, this allows us to extract from Flm for instance the 
following chiral constants 


F^ = 2Res [nLR{Q^),0] 


(76) 

(77) 


One has 


Fn=r + 


{7Z^-9)p2 

6(c2 


90k:2 {p 2 b 2 (4C(3) + 5-K^)-2{7r^- 10) P 4 ) - (-4^(3) - 5 + 

360k:6 

AfcK-2(l80C(3) + 191-41;r2) 5 - \0) (^ 4 ) 45 - \0) 


72k^ 


■-f ■ 


56 


(78) 


using Nc = 3, K = y^l.43/4GeV ci 0.6GeV and the values of the condensates (^ 4 ) = 
(-0.635 ±0.04) • 10-2 GeV^ and = (14±3) • lO^^GeV® from [4], we obtain 


- \/4099.9± 579 ±1147.8 MeV - 76 (±3)^^^ MeV , (79) 


the error in (79) are coming from the errors quoted for a and the condensates. 

In figure 4, the contributions to the final value of Fjc due to each term in (78) are reported. 
The expression for Lio is 
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Fig. 4 Relative conti'ibutions to from eq. (78). "C.T." stands for the contributions coming from the shift 
and the contact term g^p*zS{z). 


1 ^ 6Q{7i^-6)b2K^+l2Qgll52K*{-6(;{3,)-3 + 7i^) 

^ [ 5760^P _ 

r -1271:4/32^2 - 180;z:2 k- 2 (3/34 +/32Z?2 ) +360jS2fc2 (6C (3) +1) 

5760if^ 

-30ic4 {-6bl;{3) - 84^4 + {bl + %bZ)) + 540g^jS4K-2(4C(3) + 5) 

5760g2K-8 

-6g2j62 (_30(6C(3) + 1) + 15;r2 ^ ' 

5760K-8 

[7r4^3 ^2 15^2 (^3 30^4^2 + 12b(,) - ISOkt^ {b\t;{3) + 8 fo 4 ^ 2 (C(3) + 2) + 5?,be) 

^ [ 5760^P 

V,(8010C(3) +495 - 585;r2 _ 45 ^ 4 ^ 

8640;r2 

-72^(3)-12+llTT^ (^ 4 ) 5[5216C(3)+67-33;r2] 

64 if4 1792 k -6 ’ 

(80) 

then with the same numerical values used for the evaluation of F^, 

10^Lio 24-4.6-0.8 + 0.1 :ii-5.3(±l)g^t , (81) 

the error in (81) are coming from the errors quoted for a and the condensates 4. 

In figure 5, the contributions to the final value of L\q due to each term in (80) are re¬ 
ported. 

The values of Fjc and Lio are quite reasonable face with our model hypothesis, if one 
admits the usual 30% error coming from Large-Nc QCD limit, compared to the range of 
variation of in the chiral limit, 66 < F^ < 84 MeV, and compared to the value lO^Lio = 


^Let us notice that if we had made the choice for fj ~ 0.9 GeV^ as in [18] the values obtained would have 
been 80 MeV and lO^Ljo —6.2 which remain quite acceptable too. 
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Fig. 5 Relative contributions to Ljo from eq. (80). "C.T." stands for the contributions coming from the shift 
gjjSo and the contact term g^j5*zS(z)- 


—5.3 ± 0.13 according to [43] and references therein. We can also notice here that the value 
of is strongly related to the value of and then to ct = as illustrated on figure 6. 



Fig. 6 F„ variation according to a 


The relative contributions to Fj^ and Lio from the gluon and the quark condensates are 
consistent with previous evaluation [44, 45]. 


4.3 Comments on the duality violation 

The so-called duality violation phenomena [46-50] is a manifestation of the loss of the 
quark-hadron duality, i.e., the two descriptions of the same Green function of QCD, the one 
in the Minkowski region (Re ^ 0) as a collection (series) of bound states (resonances) 
and the one valid in the (deep) Euclidean region (Reg^ 0), are not related by analytic 
continuation. It implies that each expression is valid only in their own region. This two 
regions in the complex plane are disconnected and separated by the so-called Stokes’ 
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line. Until now, nobody knows where the QCD Stokes’ line is, but one suspects that it could 
he the imaginary axis. 

Our model provides analytic expressions for the vectorial and axial correlators, it is 
natural to look for such duality violations, since it is easy to recover the different behaviours 
of the functions both in Minkowski and Euclidean regions. In our analysis a crucial role is 
played by the second Weinberg sum rule of Tim, whose importance in duality violations of 
QCD sum rules has been studied, also at a quantitative level in [51]. 

As we already said in section 2.4, Tly^ corresponds to the model exposed hy the authors 
of [27, 37] and their discussion about the duality violation is still valid in our case, since 
our work is a "natural" extension including higher derivatives of the Digamma function ij/ 
in (49). 

The duality violation A = (developed in powers of 6 ) is then 




(82) 


where the functions are given in the Appendix B. is not 
since. 


is not vanishing in our model. 



(83) 


where for simplicity, a = —Therefore our model has a natural implementation for 
duality violations that will be investigated in a future work. 

5 Conclusions 

There have been several models describing QCD properties of the two point vectorial and 
axial Green function. Already an interpolation among low energy chiral properties and per¬ 
turbative QCD is good [1]. The intermediate region could be phenomenologically matched 
with a tower of resonance states equally spaced (a la Regge) [4, 6, 24, 27]. Indeed, there are 
excellent dynamical reasons that Regge trajectories are dynamically generated. Holographic 
QCD gives us a fundamental theoretical tool with the SW model to start directly from a 
theory where Regge trajectories are analytically implemented. The original model [18] had 
several difficulties: the absence of a satisfying description of chiral symmetry breaking and 
of the axial sector, wrong OPE. For some further attempts to cure these problems one can 
refer to [52, 53]. 

It is natural to modify opportunely the SW model to comply with the OPE. While it 
was done already for the HW model [45, 54], it is a novelty for the vector correlator in the 
SW model: we obtain the solution for the vector field fy profile and the analytic expression 
TI^\q^) in eq. (49) in terms of the Digamma function iff and its derivatives, from the set 
of differential equations in (34) (n is related to the order of OPE that we are addressing) - 
Appendix A and Appendix B are dedicated to a description of the calculation. As a result, 
we support previous works [24, 27, 50]. 

The coefficients of the polynomial of the dilaton profile (b^’s), fixed by the requirement 
of a correct vectorial OPE, appear also in the differential equation (63) of the axial vector 
field /a, whose expression extends the analogous ones in Ref. [10, 11, 13, 28] with the 
presence of a scalar field in the bulk with a non-trivial vacuum profile, v(z). Several issues 
have to be understood in connection also for the determination of the profile v(z): chiral 
symmetry breaking, pion pole, the correct axial spectrum and correct OPE. A coherent, 
complete and phenomenologically consistent solution emerges: due the phenomenologically 
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observed relation Sif^ = the solution for the axial in eq. (63) can be obtained from the 
vectorial one as nA{Q^) = +4if^) + "corrections". This solution generates the pion 

pole , the correct axial spectrum and axial OPE, if the profile of the vacuum, v(z), is not only 
a polynomial (jSj: in eq. (C.24)) but contains also a Dirac 5-function term, i.e. a boundary 
term for the axial field, which is needed in order to comply with the axial OPE. 

In Appendix E we show an original pictorial representation of the solutions obtained, 
both for vectors and axial vectors, up to to the inclusion of the term z® in the polynomial 
part of the dilaton. 

The necessity for simultaneous local and non-local (the Dirac 5 function contribution) 
chiral symmetry breaking mechanism has been already noticed in the literature; for instance 
in Ref [55] in the context of HW with chiral symmetry broken by boundary conditions an 
additional scalar field with canonical AdS dimensions was added to describe the pion mass. 
In Ref [52] the behaviour <P ^ of the scalar field in the UV potentially has our same 
feature: the presence of an unwanted dimension-two condensate. 

Our analytic solution for Fjc and L\q respectively in eq. (78) and (80) are very success¬ 
fully phenomenologically and show how these parameter are linked to Regge spacing and 
QCD condensates; compared to previous literature [4, 6, 24, 27] our results are fully analyt¬ 
ical. 

We think there are many applications. Eor example, our model is the fitting framework 
to analyse the duality violations since they appear naturally as the analytic continuation 
properties of the Digamma function (and its derivatives) cf. eq. (83). We enforce and extend 
here properties established from ad hoc models in Ref. [27, 37, 50]; in our case these prop¬ 
erties are based on a (5D) Lagrangian description and analytic results for the solutions of 
the corresponding equations of motion. 
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Appendix A: Integral representation of the Green function 


The solution of the unperturbed equation of motion {'Sfy^'^ = 0) could be expressed with 

( Ci') 

generalized Laguerre polynomials L„ or by an integral representation [18, 36] 




4k^ 




Au uJi? 


«=o ^ + («+ 1) 

The associated Green function is given by 


exp 


\ —u 


2 2 
K^Z 


Gv(Q^-,x,y) = — 


1 


n+\ 


4i? 


-f (n-l- 1) 


(A.1) 


(A.2) 
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such that one has explicitly 


fv\(f,x) = - limwQ{z)fv\Q^,z)d^G{Q^;x,z) , (A.3) 


z —>0 


since Li^^(O) = n + 1. 

It is more useful for our calculations to have an integral representation for the Green 
function too. One can express the Green function as an integral by using the Poisson Kernel 
[57] 


y ^>±1) 


= (xyt) '^(1 — t) * exp 


1 -r 


(x + y) 




where is the Bessel function, |t| < 1 and a > — 1, and since a = 1 in our case, writing 


1 


+ («+!) 


/•i e2 

: / dr tJi? 

Jo 


(A.5) 


then 




Gv{Q^;x,y) = -— / dtt^y -^^ 

2 Jo „to «+l 


xy n , f4i? 2 
— / dt -exp 

2 Jo l-t ^ 


1 -t 


K\^+y^) 


Ii iK'^xy 


2 .... Vi 


l-t 


(A.6) 

(A.7) 


Appendix B: Analytic expression for the vectorial two point function 

The procedure that we use for the evaluation of the analytic expression for the two point 
function is quite simple. We first use the integral representations (A.l) and (A.7), then we 
perform the integrations over the exponentials and Bessel’s functions in the integrand and 
then there remain only integrals of rational functions. As an example let focus on Fly'^ 
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calculation, 


nv\Q^) = [ dxwo{x)xdj,f^°\Q^,x)limwo{z)fv\Q^,z)Gv{Q^;z,x) 

gj Q Jo z-^o 

= fv\Q^^^) ^xfv\Q^,x) 

gi Jo 

T I / dxe^'^ ' / dudv ‘vJl? 

g 2 q2 \^4k-2; 


(B. 8 ) 

(B.9) 
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22^00 8K^g? 


/ 4if^\ 1 / 4if^\ 


V GV 3 V G^ y 


M V 

1 -^ 1 - 

1 — w 1 — V 


2 2 


(B.IO) 

(B.ll) 

(B.12) 

(B.13) 

(B.14) 


Each step can be applied to the evaluation of any higher order contributions. One always 
first performs the integrals over the hfth dimension ( [0,o°[ ) as in equation (B.IO), then 
one remains with multiple integrals between [ 0 , 1 ] of the parameters entering in the integral 
representation of and Gy - cf (A. 1 ) and (A.7)- as in (B. 11 ). For higher orders, this takes 
a general form which can be expressed as hypergeometric functions. 


[ dMi •••dM„ m "‘(1 - Ml)'’* •• •«“"(! - Ml)'’"—-i-- 

Jo (1 -Ml •••M,,)'^ 


— b(1 + ai, 1 + hi) • • - 8(1 + fl„, 1 + h„) „+iF„ 


C, 1 + Ml , . . . , 1 + 

2 + ai + fci ,...,2 + a „ + fo „ 


(B.15) 


where b is the Euler’s beta function: B(x,y) = and the hypergeometric function is 


«+iF« 


C, 1 + Ml, . . . , 1 + 

2 + Mi + &!,...,2 + Mn+hn 


j \ _ {^)kO d~ 44l}k‘' ' 0 44n)k 

' ^oO)k{2 + ai+bi)i,---{2 + a„+b„)i, 


(B.16) 

where the Pochhammer symbols are dehned as {a)k = r{a + k)/r{a). The hypergeometric 
functions ^F^ of argument 1 have the particularity that their coefficients are function of ^ ± 
1, so that in this case the series in (B.16) can be summed in terms of Digamma functions 
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and rational function [58] as summarized in (49), viz. 
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Appendix C: Analytic expression for the axial two point function 


We shall proceed in analogy with what we did in section 2.4. We expand, 

fA = /f’ + e/i'^ + eVf + eVf + ^(e") 

and the associated axial two-point function, 

Ha = + ff { e ^). 


(C.21) 
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The expression (71) is also expanded in powers of d. 




= pQ+p2z^e+P4z^e^+p*z5{z). 


(C.23) 


(The absence of powers of 9 for the /3* term is due to its invariance under the rescaling 

z ^ VOz.) 

The equation of motion for the axial field (63) becomes 
fA + d~^r\w{z)\dyfA — [Q^ + g\fk))fA = ^5 [Piz^d + I34z‘^6^ + p* z5{z )]/a , (C.24) 

with the boundary conditions 

/W(e2,0) = 5o.„ and (C.25) 

The zero-th order solution of (C.24), neglecting the contribution from the jS* term, is 
given by (simply the implementation of the shift on the integral representation (46) ) 
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and the associated Green function. 
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The P* term contributes only to the 9^ term, and as it was the case for the vector correlator 
in Sect. 2.4, also for the axial one, we obtain, order by order in 9 : 

<?5 *^0 

X E IglPiin-kA’'-’'^-^ djf\Q\x)] . (C.28) 


Using the same methodology one obtains the following contributions proportional to the 
Pi, coefficients in (C.28) for the axial two point functions, 
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Appendix D: Corrections to the mass spectrum 


In this Appendix, we want to show explicitly our assumption that Regge trajectories are kept 
in our approach and modified order by order in d by sub-leading corrections in n. Indeed, the 
parameter 9 introduced in Sect. 2.3 acquires a physical meaning as a true perturbative pa¬ 
rameter, when we try to quantify corrections to the leading Regge mass spectra. We develop 
the considerations we did at the end of Sect. 2.6 for the vector resonance mass spectrum. 
The same argument would hold for the axial resonances as described in Sect. 4. 

As we showed in Sect. 2.4, the vector two-point function nv{Q^) receives corrections 
in 9 eq. (51) as 

nviQ^) = 'tnl!‘\Q^)9’^, (D.31) 

k=0 

which generically take the form 




h\h\ (e^+4K:2«)^ ’ 


(D.32) 


showing the appearance of multiple poles of increasing order the greater the powers of 9 
considered in eq. (49). 

The question arises whether it is possible to rewrite (D.32) in a way that preserves the 
Regge behaviour, maybe with the cost of modifying the value of the slope parameter. In 
fact, it is possible to answer in the affirmative, provided the value of 9 is small enough. It is 
possible to show that (D.32) can be generically rewritten such as 


nviQ^) 


y Fv{n,9)^ 
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(D.33) 


where 

a{n, 9) = 4ic^ + _|_... 
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(D.34) 

















27 


corresponding to corrections to the original Regge slope, a = 4 k^, and to the residue. The 
equivalence between the two expressions (D.32) and (D.33) is non trivial, and relies on a 
number of algebraic relations satisfied by the coefficients Af'^ (n) that reduce the number 
of independent parameters thus allowing the equivalent representation (D.33). Thus, for a 
given resonance term, i.e. a fixed value of «, and a given value of the power of 9 considered, 
the departure from the original Regge slope can be made parametrically small by choosing 
a correspondingly small enough value of 9. Notice that coefficients of higher power in 9 are 
expected to be sizeable or even greater than the lower ones, following, at most, an asymptotic 
series behaviour. 

In this context, the knowledge of the complete OPE series would drive to a complete 
resummation of the Regge slope since for the vector correlator there is a one by one corre¬ 
spondence between the dilaton modification coefficients and the OPE coefficients eq. (56). 
We have illustrated here that the partial information of few OPE terms can be treated in our 
model such that the QCD properties on the Minkowski region are still satisfied, in other 
words, we assume Regge trajectories even if it is realized dynamically and impose a correct 
OPE. 


Appendix E: Diagrammatic representations of the results 

In order to simplify the understanding of the calculation procedure, we can represent graph¬ 
ically the result by constructing diagrammatic rules. Let us define the bulk-connectors and 
boundaries-connectors, as illustrated on tables 1 and 2. One can associate to each connector 
an integral and just by multiplying each of them, we obtain the corresponding contribu¬ 
tion. The last rule is that to obtain all the contributions to TIy'\ one needs to collect all the 
possible ways to have a path with n legs starting from the boundary-connector. 

One can see all the possible combinations on figure 7 for the vectorial case and on figure 
8 for the axial case. For example, it is easy to recover eq. (B.8) using this rules. 



F{x) 
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k times 



= l’2k f dxwa{x)Gv{Q^;z,x)J‘d^fl^°\Q^,x) 
Jo 

= b2k f dxwo{x)Gv{Q^-,Z,x)x!‘d:cF{x) 

Jo 

= g5& ^ 4xwo{x)GA{Q^-,z,x)x'‘jf\Q^,x) 
= g\^k f dxwo(x)GA{(f;z,x)x'‘F{x) 

Jo 


Table 1 Diagrammatic rules of construction: the bulk-connectors. On the right one has the corresponding 
operator. 
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= ^“2 / dzwo(z)4°’(e^,z)5(z)aF(z) 

§5 •'0 

= “ “2 ^“'9, ^0 (2^ z) SjF (z) 

= ”“2 f dzwo(z)/i“'(e^,z)5(z)5zF(z) 
§5 20 


Table 2 Diagrammatic rules of construction: the boundary-connectors. On the right one has the correspond¬ 
ing operator. 



Fig. 7 All graphical contributions to the vectorial two-point function up to 
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from the shift are inside the dotted blue rectangle {i.e. the path containing only blue lines). 


























30 


References 

1. S. Peris, M. Perrottet and E. de Rafael, "Matching long and short distances in large N(c) 
QCD", JHEP 9805, 011 (1998) 

2. P.D.B. Collins, "Regge theory and particle physics", Phys. Reports 1 (1971) 103. 

3. G. ’t Hooft, Nucl. Phys. B75 (1974)461; 

C.G. Callan et al., Phys. Rev. D13 (1976) 1649; 

M.B. Einhorn, Phys. Rev. D14 (1976) 3451. 

4. M. Golterman and S. Peris, "Large N(c) QCD meets Regge theory: The Example of 
spin one two point functions" JHEP 0101, 028. 

5. A. B. Kaidalov, "Regge poles in QCD", In Shifman, M. (ed.): At the frontier of particle 
physics, vol. 1 603-636 

6. M. Shifman and A. Vainshtein, "Highly Excited Mesons, Linear Regge Trajectories and 
the Pattern of the Chiral Symmetry Realization", Phys. Rev. D 77, 034002 (2008) 

7. J. M. Maldacena, "The large N limit of superconformal field theories and supergravity", 
Adv. Theor. Math. Phys. 2, 231 (1998). 

8. S. S. Gubser, I. R. Klebanov and A. M. Polyakov, "Gauge theory correlators from non- 
critical string theory", Phys. Lett. B 428, 105 (1998) 

9. E. Witten, "Anti-de Sitter space and holography", Adv. Theor. Math. Phys. 2, 253 (1998) 

10. J. Erlich, E. Katz, D. T. Son and M. A. Stephanov, "QCD and a holographic model of 
hadrons", Phys. Rev. Lett. 95, 261602 (2005) 

11. L. Da Rold and A. Pomarol, "Chiral symmetry breaking from five dimensional spaces", 
Nucl. Phys. B 721, 79 (2005) 

12. J. Him and V. Sanz, "Interpolating between low and high energy QCD via a 5-D Yang- 
Mills model", JHEP 0512, 030 (2005) 

13. D. T. Son and M. A. Stephanov, "QCD and dimensional deconstruction", Phys. Rev. D 
69, 065020 (2004) 

14. T. Sakai and S. Sugimoto, "Low energy hadron physics in holographic QCD", Prog. 
Theor. Phys. 113, 843 (2005) 

[15] 

15. S. J. Brodsky, G. F. de Teramond, H. G. Dosch and J. Erlich, "Light-Front Holographic 
QCD and Emerging Confinement", arXiv: 1407.8131 [hep-ph]. 

16. A. A. Migdal, "Series Expansion for Mesonic Masses in Multicolor QCD", Annals 
Phys. 110, 46 (1978). 

17. M. Shifman, "Highly excited hadrons in QCD and beyond", hep-ph/0507246. 

18. A. Karch, E. Katz, D. T. Son and M. A. Stephanov, "Linear confinement and 
AdS/QCD", Phys. Rev. D 74, 015005 (2006). 

19. O. Andreev, "l/q**2 corrections and gauge/string duality", Phys. Rev. D 73, 107901 
(2006) 

20. O. Cata, "Towards understanding Regge trajectories in holographic QCD", Phys. Rev. 
D 75, 106004 (2007) 

21. F. Zuo and T. Huang, "Comments on the two-dimensional power correction in the soft 
wall model". Chin. Phys. Lett. 25, 3601 (2008) 

22. C. Csaki and M. Reece, "Toward a systematic holographic QCD: A Braneless ap¬ 
proach", JHEP 0705, 062 (2007) 

23. G. Veneziano, "Constmction of a crossing - symmetric, Regge behaved amplitude for 
linearly rising trajectories", Nuovo Cim. A 57, 190 (1968). 

24. M. Golterman, S. Peris, B. Phily and E. De Rafael, "Testing an approximation to large 
N(c) QCD with a toy model", JHEP 0201, 024 (2002) 



31 


25. C. A. Dominguez, "Pion form-factor in large N(c) QCD", Phys. Lett. B 512, 331 (2001) 

26. S. S. Afonin, A. A. Andrianov, V. A. Andrianov and D. Espriu, "Matching Regge theory 
to the OPE", JHEP 0404, 039 (2004) 

27. O. Cata, M. Golterman and S. Peris, "Duality violations and spectral sum rules", JHEP 
0508, 076 (2005). 

28. H. J. Kwee and R. F. Lebed, "Pion Form Factor in Improved Holographic QCD Back¬ 
grounds", Phys. Rev. D 77, 115007 (2008) 

29. T. Gherghetta, J. I. Kapusta and T. M. Kelley, "Chiral symmetry breaking in the soft- 
wall AdS/QCD model", Phys. Rev. D 79, 076003 (2009) 

30. M. Beneke, "Renormalons", Phys. Rept. 317, 1 (1999) 

31. G. ’t Hooft, "A Planar Diagram Theory for Strong Interactions", Nucl. Phys. B72 (1974) 
461. 

32. E. Witten, "Baryons in the 1/N Expansion", Nucl. Phys. B160 (1979) 57. 

33. J. Beringer et al. [Particle Data Group Collaboration], "Review of Particle Physics 
(RPP)", Phys. Rev. D 86, 010001 (2012). 

34. P. Masjuan, E. Ruiz Arriola and W. Broniowski, "Systematics of radial and angular- 
momentum Regge trajectories of light non-strange qq-slates,", Phys. Rev. D 85, 094006 
( 2012 ) 

35. M. Shifman, A. Vainshtein and V. Zakharov, "QCD and Resonance Physics. The rho- 
omega Mixing", Nucl. Phys. B147 (1979) 385; 447. 

36. H. R. Grigoryan and A. V. Radyushkin, "Structure of vector mesons in holographic 
model with linear confinement" Phys. Rev. D 76, 095007 (2007). 

37. M. Jamin, "What two models may teach us about duality violations in QCD", JHEP 
1109, 141 (2011). 

38. E. de Rafael, "Large Nc QCD and Harmonic Sums" Pramana-joumal of physics, Indian 
Academy of Sciences, Vol. 78, N6 June 2012 pp. 927-946. 

39. J. Mondejar and A. Pineda, "1/N(c) and 1/n preasymptotic corrections to Current- 
Current correlators", JHEP 0806, 039 (2008) 

40. L. Cappiello, G. D’Ambrosio and D. Greynat, "The pion wave function in a Soft Wall 
Model", To be published 

41. P. Cox and T. Gherghetta, "A Soft-Wall Dilaton", JHEP 1502 (2015) 006 

42. P. Colangelo, F. De Fazio, J. J. Sanz-Cillero, F. Giannuzzi and S. Nicotri, "Anomalous 
Ay*V vertex function in the soft-wall holographic model of QCD", Phys. Rev. D 85, 
035013 (2012) 

43. G. Ecker, "Facets of chiral perturbation theory", Nucl. Phys. Proc. Suppl. 245 (2013) 1 

44. E. de Rafael, "Chiral Lagrangians and kaon CP violation". In *Boulder 1994, Proceed¬ 
ings, CP violation and the limits of the standard model* 15-84, and Marseille Cent. 
Theor. Phys. - 95-P-3161 (95/01,rec.Feb.) 71 p. (507515) 

45. J. Him, N. Rius and V. Sanz, "Geometric approach to condensates in holographic QCD", 
Phys. Rev. D 73, 085005 (2006) 

46. M. A. Shifman, Theory of pre-asymptotic effects in weak inclusive decays. In Proc. of 
the Workshop Continuous advances in QCD, ed. A. Smilga (World Scientific, Singa¬ 
pore, 1994). 

47. M. A. Shifman, Recent progress in the heavy quark theory. In Particles, Strings and 
Cosmology, eds. J. Bagger et al. (World Scientific, Singapore, 1996). 

48. B. Chibisov, R. D. Dikeman, M. A. Shifman, and N. Uraltsev, Operator product ex¬ 
pansion, heavy quarks, QCD duality and its violations, Int. J. Mod. Phys. A12 (1997) 
2075-2133,. 




32 


49. B. Blok, M. A. Shifman, and D.-X. Zhang, An Illustrative example of how quark hadron 
duality might work, Phys. Rev. D57 (1998) 2691-2700. 

50. M. A. Shifman, Quark-hadron duality. Published in the Boris Ioffe Festschrift ’At the 
Frontier of Particle Physics Handbook of QCD’, ed. M. Shifman (World Scientific, 
Singapore, 2001). 

51. M. Gonzalez-Alonso, A. Pich and J. Prades, "Pinched weights and Duality Violation in 
QCD Sum Rules: a critical analysis", Phys. Rev. D 82, 014019 (2010) 

52. B. Batell and T. Gherghetta, "Dynamical Soft-Wall AdS/QCD", Phys. Rev. D 78, 
026002 (2008) 

53. G. F. de Teramond and S. J. Brodsky, "Light-Front Holography and Gauge/Gravity Du¬ 
ality: The Light Meson and Baryon Spectra", Nucl. Phys. Proc. Suppl. 199, 89 (2010) 

54. L. Cappiello and G. D’Ambrosio, "On the Evaluation of Gluon Condensate Effects in 
the Holographic Approach to QCD", Eur. Phys. J. C 69, 315 (2010) 

55. O. Domenech, G. Panico and A. Wulzer, "Massive Pions, Anomalies and Baryons in 
Holographic QCD", Nucl. Phys. A 853, 97 (2011) 

56. P. Gambino, "B semileptonic moments at NNLO", JHEP 1109, 055 (2011) 

57. E. D. Rainville, "Special Functions", Macmillan, NewYork, 1960. 

58. J. E. Gottschalk and E. N. Maslen, "Reduction formulae for generalised hypergeometric 
functions of one variable", J. Phys. A. Math. Gen. 21 (1988) 1983-1998. 




